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We present a simple one-dimensional model with molecular 
interactions favouring the formation of clusters with a defined 
optimal size. Increasing the density, at low temperature, the 
system goes from a nearly ideal gas of independent molecules 
to a system with most of the molecules in optimal clusters, in 
a way that resembles the formation of micelles in a dilution of 
amphiphilic molecules, at the critical micellar concentration. 
Our model is simple enough to have an exact solution, but it 
contains some basic features of more realistic descriptions of 
amphiphilic systems: molecular excluded volume and molec- 
ular attractions which are saturated at the optimal cluster. 
The comparison between the exact results and the mean field 
density functional approximation suggests new approaches to 
study the more complex and realistic models of micelle forma- 
tion; in particular it addresses the long-standing controversy 
surrounding the separation of internal degrees of freedom in 
the formulation of cluster association phenomena. 



I. INTRODUCTION 

The study of systems with strong effects of molecular 
aggregation sets problems of fundamental and applied 
interest. One of the most important examples are the 
dilute solutions of amphiphiles in waterEriJ. The polar 
head in the molecule attracts the water, which is re- 
pelled by the hydrocarbon tails in the same molecule, 
and the frustration of these opposite tendencies leads 
to the formation of micelles, droplets with hundreds or 
thousands of molecules, which have their polar heads on 
the surface (in contact with the water) and the interior 
filled with the hydrocarbon tails. The radius of the mi- 
celles is limited by the length of the molecular chains, 
so that the droplets cannot grow far beyond it as spher- 
ical drops. From a theoretical point of view, micellar 
solutions could be regarded as extreme cases of non-ideal 
solutions, with strong molecular correlation described as 
microscopic clusters (micelles). However, the large num- 
ber of molecules in each micelle makes impracticable t 
use of the virial expansion as with ordinary dilution; 
The alternative is to split the description of the systems 
into two levels: first the description within some approx- 
imate scheme, of a single cluster of k molecules to get its 
internal grand potential energy, fife , and then the whole 
system is regarded as an ideal mixture, with the density 
of clusters proportional to the internal Boltzmann factor 



with (3 = (fcf,T)~^. There remains controversy about the 
proportionality factor, which should provide the units 
of inverse volume to in Eq. (^, and how liis fac- 
tor may depend on the cluster size and masaJ'El; see 
also the series of articles by Reiss et alu, and references 
therein. The indeterminacy in this factor is directly as- 
sociated to the problcan of how to define the internal 
grand-potential energyS; when the process of micelle for- 
mation is described as a chemical reaction the problem 
is transferred to thft-psual indeterminacy in the standard 
chemical potentialaa. The problem may be avoided as- 
suming that the dependence of pk with the temperature 
and the total concentration of amphiphiles is dominated 
by the exponential factor. For fi^ > micelles of size k 
are scarce, while for 17^ < they should be abundant; 
the critical micellar concentration (CMC) is then asso- 
ciated with conditions at which the minimum value of 
il/c, for the optimal size, n, becomes zero. This semi- 
quantitative analysis provides useful information from 
empirical thermodynamic approximations for Ofc, with 
surface and volume contributions, which may include the 
qualitative effects of molecular shaped and the role of the 
flexible chainsotl. However, attempts to develop a quan- 
titative link between approximated (mean field) descrip- 
tions of an aggregate at microscopic le«gl, and the density 
Pk of these aggregates at equilibriumoB, requires consis- 
tent microscopic approximations for the prefactor in Eq. 
(|l|). In this paper we explore the problem with a sim- 
ple model for molecular aggregation that may be solved 
both exactly and with a mean field density functional 
approximation XMFA) comparable to that used in more 
realistic modelaS. The comparison allows us to establish 
a workable approach linking the two levels of statistical 
description in these systems, to be used in more sophis- 
ticated models which do not have exact solution. 



II. SIMPLE MODEL FOR MOLECULAR 
AGGREGATION 

Our model for molecular aggregates is a system of hard 
rods, with length u, moving along the X axis. Consecu- 
tive rods have positions x^+i > Xi + a, and each rod has 
a discrete internal degree of freedom ^ = 1, n. Besides 
the hard-core repulsion there is an interaction energy be- 
tween nearest neighbour pairs, 4>{xi+i — Xi,^i+i — ^i) re- 
stricted to the range a < x^+i — Xi < 2a. As a simple 
model with this behaviour we take 
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{x — a){x — 2a) 
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if ^ = — S.i = ±1 and zero otherwise. The dimen- 
sionless parameter K gives the inverse temperature, in 
units of the interaction strength. The optimal aggregate, 
represented in Fig. ||, is made of n rods with consecutive 
values of from = 1 to ^„ = n, keeping relative dis- 
tances between neighbours close to Iq = 3tT/2, to be at 
the minimum of (f>{x,^). The non-cyclic character of the 
internal variable ^ precludes the growth of larger aggre- 
gates, playing (in a rather different way) the same role as 
the chain-packing constrains in real amphiphiles. Small 
fluctuations around the optimal distance x^+i — Xi = Iq 
still have a cohesive energy and provide some internal 
entropy for the aggregate. The aggregates compete with 
fully disordered configurations, which maximize entropy 
but do not get cohesive energy from Eq. (S). 




FIG. 1. Sketch of an optimal molecular aggregate in our 
one dimension model with n — 4. The hard rods are separated 
by the optimal distance, Iq, and the internal variable takes 
consecutive values from = 1 to — n as indicated by 
the numbers in each rod in (a). In the lower representation 
(b), the internal variable is represented by an angular variable 
with half-clock n possible values, to give a more intuitive idea 
of how the non-cyclic character of prevents the grow of 
aggregates with more than n rods. 



A. Exact results 

The exact properties of the model are easily obtainedS 
in a statistical ensemble with constant temperature, T, 
and pressure, p, in terms of a (n x n) transfer matrix 
with elements, Mj,k{x) = exp[—l3{(j){x, k — j) — px)] for 
j,k — l,..,n. The integral of Mj_k{x) from a; = 1 to 
infinity gives the elements of the transfer matrix for the 
internal variable, Qj,fc, and its largest eigenvalue Xmax 
provides the thermodynamic properties of the system, 
through the chemical potential l3^{T,p) = — ln(Amaa;), 
and the density of rods, p{T,p) = {dp/dp)T- The 
isotherms K = 2Q for n = 2 and n = 12 presented in 



Fig. II A, give (3iJ. versus pa/n (with logarithmic scale). 
Very diluted systems present ideal gas behaviour, with 
curves of slope 1 in the figure. When the molecular as- 
sociation becomes important the curves change to have 



slope 1/n, which corresponds to an ideal mixture of n- 
molecule clusters, until the packing effects show up and 
the slopes increase. 




logi„(pa/n) 

FIG. 2. Isothermal representation of the chemical poten- 
tial versus the molecular density in our model for if = 20. 
The full lines are the exact results for n = 2 and n = 12. Bro- 
ken lines are the results with the mean field approximation. 
The dotted line is the ideal gas limit. 



Any property measuring the correlation structure can 
also be obtained from M{x) and Q. In particular, the 
density of aggregates, which is a property of n-rod cor- 
relations, is 
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where the matrix W depends on our exact definition for 
the aggregate. In the following we define the full ag- 
gregate (or micelle in a rather loose sense of the word) 
as any group of n rods with a < Xi+i — xi < 2(j, and 
^i+i — = 1, starting with ^ = 1, so that every bond 
gets some cohesion energy from Eq. (|). With this defi- 
nition we have 



= (5j,i(5fc,„ / dxMi.i+i{x). 
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Slightly different definitions of micelle would lead to other 
forms for W, but they would change p„ very little, as long 
as they include the optimal and nearby configurations. 

In a similar way we may calculate the density, pk , of in- 
complete aggregates^ with k < n rods, and the density pi 
of isolated r ods, which do not get any cohesion from Eq. 
(^. In Fig. [I A we present k pk/ p, i.e., the proportion of 
rods forming part of cohesive clusters of different sizes, as 
a function of the total density p, along an isotherm. In 
qualitative agreement with real systems of amphiphilic 
molecules, there is a low concentration regime, with very 
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few micelles, where p ~ pi, and a higher concentration 
regime, dominated by clusters of the preferential size n, 
so that p ~ n p„ . The density at the transition between 
the two regimes represents the critical micellar concen- 
tration (CMC, again in a loose sense of the term), which 
may be defined at the point where pi = n p„, so that 
there are as many isolated molecules as molecules form- 
ing part of complete micelles. In the neighbourhood of 
the CMC there are significative amounts of incomplete 
aggregates, with size between 2 and n — 1. The CMC 
correspo nds to the region where the slopes of the curves 
in Fig. [I A change from 1 to 1/n. These changes are 
smoother in our simple, one-dimensional model than in 
real amphiphilic systems, but as the CMC is never a true 
phase transition, the use of a one-dimensional model does 
not change the quantitative behaviour of the system. 
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FIG. 3. Proportion of rods forming part of aggregates of 
different size in our model for n — 5, as a function of the 
molecular density along the isotherm K — 60. The curve 
with k = 1 corresponds to isolated molecules, the curve with 
k = n gives the proportion of rods in full aggregates. The 
critical concentration region is located around the point where 
pi = n p„, and has high relative concentration of incomplete 
aggregates, with 2 < fc < n — 1. 



At low temperatures {K ^ 1) the asymptotic analysis 
of our exact results gives the CMC at 
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independent of the optimal cluster size n. In Fig. II A| we 
compare Eq. (^) with the exact results obtained from the 
numerical solution of Eq. (|[ ^). The deviations from the 
asymptotic limit at low K (i.e., high temperature) repre- 
sent the effects of the molecular packing, which force the 
relative distances between neighbour rods below the opti- 
mal value Zqi with a compromise between the tendency to 
build the optimal aggregates and the packing constrains 
(which push the rods towards the distance a given by 



the hard core). In real systems of amphiphilic molecules, 
this compromise leads to very complex naiiase diagrams, 
with lamellar, hexagonal, cubic phasesEJ'U. In our one- 
dimensional system with short range interactions, how- 
ever, the exact result cannot include any phase transi- 
tions. 




FIG. 4. The critical concentration in our model repre- 
sented in terms of the parameter K {or adimensional inverse 
temperature) versus the chemical potential /3/i. The full lines 
are the exact results for n = 2 and n = 12. They share the 
same asymptotic limit at low temperature [K 1), as given 
by Eq. (^ and represented by the dotted line. The results 
with the mean field approximation are given by the dashed 
(n=12) and the dot-dashed (n=2) lines. They have an asymp- 
totic shift respect to the exact results which decreases with 
the aggregate size n. 



B. Mean field approximation 

We now consider the same model with a density func- 
tional approximation, to describe the aggregates as meso- 
scopic self-maintained inhomogeneities in a homogeneous 
bulk phase, corresponding to local minima of the grand- 
potential energy as a functional of the distribution of rods 
with each value of the internal parameter, p(x,^), 

^[pi^^O] = ^hr[p{x)] -pi dxp{x) + 
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dxdx' 4>{x' — x,k — i)p{x, j)p{x' , k), (6) 



x<x' 



were T^r ^ the exact free energy density functional for 
hard rodstil, which depends only on p{x) = 
and the second line is the mean field (mf) approximation 
for the interaction Eq. (||). This treatment of the model 
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is equivalent to what can be done for more realistic mod- 
els in three dimensionaS. As in those cases, n[p{x,S,)] 
in Eq. always has a local minimum, with value fio, 
for a homogeneous density distribution p(a:,^) = p/n, 
in which, for our model, the contribution of </>(a:, ^) van- 
ishes, and the system follows the thermodynamics of a 
hard rod system. However, at low temperature other lo- 
cal minima may appear, with value Qo + ^k, in which the 
density distribution p(x,£,) is made of a series of k n arrow 
peaks, over the uniform background density (Fig. p^I B| ). 
The distance between the peaks and the dependence with 
^ corresponds to configurations close to the optimal ag- 
gregate. This mean field representation is equivalent to 
what is found in more realistic models of micelles in three 
dimensions. As in those cases, the value of r2[/9(a;,^)] is 
degenerate with respect to a rigid translation of the mi- 
celle, which may be located at any point over the uni- 
form background. The density of this background rep- 
resents within the mean field approximation the density 
of monomers, pi, and it may be used as a way to con- 
trol the chemical potential, p, = phriPi^T). This mean 
field approximation neglects the entropy associated with 
the position of the micelle, and the second step of the 
statistical analysis corresponds precisely to include it to 
obtain the aggregate density, pk, from their grand po- 
tential excess flk, like in Eq. (|^). Our simple model, 
with exact results for pk, allows us to explore workable 
approximations for the proportionality factor in Eq. (|^). 
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FIG. 5. The mean field representation of an aggregate as 
an inhomogeneous density distribution which is a local mini- 
mum of the grand potential energy in our model with n — 2, 
for K — 16 and bulk density pi = 0.02cr~^. The full line is 
the total density with two narrow peaks, and the dashed line 
is the internal order parameter ^{x) which is nearly saturated 
at the value ^ = 1 on the left peak and ^ = 2 on the right 
one. 

The degeneracy of il[p{x,£,)] with respect to the po- 
sition of the micelle also implies the existence of local 



minima with two or more micelles well separated along 
the system, and with an excess grand potential which is 
ilk times the number of micelles. Some weak effective at- 
traction may appear between micelles, mediated by the 
monomer density, pi, and as the chemical potential in- 
creases, the absolute minimum otfl[p{x, k)] goes to dense 
periodic structures. 

III. IDEAL MIXTURES OF MOLECULAR 
AGGREGATES 

A. General approach 

In this section we consider the generic problem of how 
to obtain the global behaviour of a system from a mean 
field description of the molecular aggregates in the limit 
of low total density, p. We use the notation appropriate 
for systems in three dimensions to get general expres- 
sions, which we later apply to our simple model in one 
dimension. We assume that there may be several types of 
molecular aggregates, and we label them by their number 
of molecules k. There is always some arbitrariness, at this 
level of description, in what is called an n-molecule ag- 
gregate, but the existence of micelles with a well defined 
size and structure makes it easy to find natural choices, 
like the one used in our simple model leading to (2-3). 
As discussed above, the use of reasonable variants would 
have very little effect in the results. 

In the limit of low total density the system may be 
considered as an ideal gas of aggregates. The partition 
function of an ideal mixture with Nk aggregates of size 
k, and a total number of molecules N = 'Ylk^-^k, is 
factorized into the contributions of each aggregate, 

k 

where Zk{T, V) is the canonical partition function of the 
fc-molecules aggregate, as a function of the temperature 
and the total volume of the system, V. This partition 
function may be factorized into a translational part, pro- 
portional to the volume, and the internal part, Z™*(T) 
which is independent of V. We make it explicit as 

Zk{T, V) = e-P^" ^^^^^^ - ^ Zl-\Tl (8) 

where we have separated the Boltzmann factor, with the 
mean energy of the aggregate Ek, and the phase space 
volume with a factor V . The k — 1 factors (for which 
Eq. (^ may be considered a definition) include the phase 
space volume associated with the relative positions of all 
the molecules in the aggregate and any rotational or in- 
ternal degree of freedom of the molecules, as well as any 
combinatorial factor. The total configuration integral 
over positions has dimensions of volume to the power 
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k, which are cancelled by the factor F'^ in the denomina- 
tor, where T is the unit cell volume in the configuration 
phase space. The matching between classical and quan- 
tum statistics may be used to fix F, e.g., for a monatomic 
ideal gas with atomic mass to, the integral over the mo- 
menta give F = A"^, in terms of the thermal wavelength, 
A = h{2nf3/m,y/^. However, within the range of valid- 
ity for the classical statistical mechanics all the results 
related to the positions of the molecules, like pk or any 
molecular correlation structure, will be invariant if F is 
multiplied by an arbitrary factor (or if the molecular mass 
TO is changed). Thus, we keep the term F*"' in Eq. (H), to 
get a dimensionless partition function Zk, but we know 
that F and m have to disappear in any final expression 
for the density of aggregates. By taking F as an arbitrary 
constant (independent of T) we neglect the molecular ki- 
netic energy, which would add (for each classical degree 
of freedom) a contribution (2/3) ~^ to the internal energy, 
without changing any relevant result here. 

The mean number of fc-molecule aggregates in the sys- 
tem, TVfe, is easily obtained from Eq. (uu in the grand- 
canonical ensemble, = —d{Zk exp{kPiJ.))/dl3Ef;, and 
the factor V in Eq. (||) provides a well defined density of 
clusters 



Pfc = ^ = — exp [-P{Ek - kfi) - fcln(F/wfe) 



The total density is then given by 
N 

fc>i 



V 



pi+^ k Pk, 



(9) 



(10) 



where we have separated the density pi of isolated 
molecules, which are not considered to form part of any 
aggregate. For a given total density, p, the chemical po- 
tential may be written as 



^ = + Am = In(pF) + Ap{p, T), 



(11) 



where we have separated the ideal gas contribution from 
the contribution of the molecular interactions. The for- 
mer takes all the dependence on the phase space unit 
F, through the form (3pid = In(pF). Substitution of Eq. 
( pi] ) into Eq. (|^) shows the perfect cancellation of the F 
factors in the density of aggregates, pk{p,T) when it is 
written as a function of the total density and the tem- 
perature. 

Consider now the description of the same system with 
a generic mean field density functional approximation for 
the grand potential energy fl[p]. There is always a rela- 
tive minimum of f2 for a homogeneous bulk density dis- 
tribution, which we associate with the density pi (T, p) 
of unaggregated molecules in Eq. (p^. There may be 
other functional local minima, for inhoniogcneous distri- 
bution functions, and each one is associated with a type 
of molecular aggregate. Compared with the uniform bulk 
phase, this aggregate has an excess of k molecules (which 
again we use to label the type of aggregate) and an excess 



Helmholtz free energy, . In the limit of a very dilute 
bulk we may represent the mean field approximation for 
the partition function on the k molecules as 



,m/-i 




(12) 



where E^''^ is the energy given by the mean field approx- 
imation, and , defined through Eq. (p^), gives the 
volume in the configuration phase space per molecule in 
the same approximation. The qualitative difference be- 
tween Eq. ( |l^ ) and the exact result Eq. (H) is the lack 
of a volume factor V. This factor arises from the overall 
position of the aggregate, but in our mean field descrip- 
tion we are neglecting the collective mode associated to 
fc-molecules correlations. The translational invariance is 
still reflected in the degeneracy of the density functional 
minimum, but it does not contribute to FjJ^^ , which is 
independent of the total volume in the system. 

We now make explicit the use of the mean field approx- 
imation for the statistics of the aggregates by using the 
exact expression Eq. (||) for pk, with the results of the 
mean field approximation for E^ and v^', the expression 
in the exponential becomes 



/3(S™-^ 



kp) - k HT/vk) = -/3F^"^ + fc/3/i (13) 



which is precisely —(30,^, , as in empirical treatments, 
but we now get an explicit approximation for the prefac- 
tor: 



Pk 



1 

7' 



where may be obtained through Eq. ( |12[ ) 



(14) 



from the 



difference between the internal energy and the Helmholtz 
free energy in the same mean field approximation. 



1 

Pi 



exp 



E 



^-/3Am(pi) 



(15) 



in terms of the bulk density of unaggregated molecules, 
the internal energy E™^ and the grand potential excess 



Within this approach we may obtain the thermody- 
namics and the structure of a system with a strong de- 
gree of molecular aggregation from a mean field descrip- 
tion of each separated cluster. For each local minimum 
of r2[p] we get a type of cluster, which at a given tem- 
perature and chemical potential, may coexist with a bulk 
of density pi{T, p) of uncorrelated molecules. The total 
density of the system is given by Eq. (p^), adding to pi 
the contribution Eq. ( p^ of each type of aggregate (i.e., 
each local minimum of fJ). The proposed expresion Eq. 
( p^ for the prefactor in Eq. (|lj), is based on the for- 
mal comparison of the exact and the mean field partition 
functions. The factor V of the free translational mode. 
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which is missing in the MFA, is added in substitution of 
an average mode with contribution t'™^. There are dif- 
ferent recipes to get Vk, and most of them are based on 
the phase space unit T (taken in terms of the thermal 
wavelength), which depends on the molecular massB'S, 
inconsistent with the use of classical statistics. An em- 
pirical recipe proposed by usO, was based on the micelle 
compressibility as a way to estimate the contribution of 
the average mode. That approximation was used in the 
description of micelles, in a three-dimensional model of 
amphiphilic systems, with sensible results. However, it 
was too sensitive to the structure of the aggregates, and 
it was difficult to use for crystal micelles with strong in- 
ternal structure. 

Our present proposal is a consistent thermodynamic 
evaluation of from the entropy of the aggregate, as 
given in the same MFA as used for the grand potential 
energy f^™'^. This requires getting the grand potential 
and the internal energy from the same approxima- 
tion. Obviously, this level of description is always based 
on an approximation; there is not an exact value for the 
prefactor in Eq. (|l|) because the exact description of the 
system (as we have done with the simple one-dimensional 
model) would not give any local minima of the grand po- 
tential energy for inhomogeneous density distributions, 
and there is not an exact value for the exponent Qn- In 
the next subsection we use our simple model of molecular 
aggregation to check the precision of this approach by di- 
rect comparison with the exact results, which is obviously 
not possible in more realistic models. 



B. Application to the one-dimensional model 

It is straightforward to use Eqs. (|l^, |l| and |l|) for 
our simple model in one dimension. The mean field ap- 
proximation Eq. (|^) has explicit separation between the 
internal energy and the entropic contributions to 
At low temperatures (K 3> 1), we may use a Gaussian 
parametrization of the density distribution: 

rn 

p{x, C) = ^ + V exp[-a,(x - x,)2], (16) 
n ^-^ 

i=l 

with Ci(f), a,; and free variational parameters 

to minimize ^[p]- The homogenous diluted bulk with 
density pi corresponds to have Ci{^) = for all i 
and With the same value of the chemical potential, 
/3/Lt(T, pi) — ln{piT/n), we may find other local minima 
with the structure of complete or incomplete aggregates, 
with k < n Gaussian peaks with the correct sequence 
of the internal variable, Ci(^) ~ (a^/Tr)^/^, if f = i and 
Ci(^) ~ otherwise. The distance between Gaussians is 
given by the minimum of the interaction potential Eq. 
(g), Xi+i — Xi = Iq = 3(t/2, and the Gaussian widths are 
ai ~ AK/a^ for the inner Gaussians and ai ~ 2K/a^ for 
those at the ends of the aggregates. 



Within this parametrization, with normalized Gaus- 
sian peaks, the excess of internal and Helmholtz free en- 
ergies of an aggregate of size k are obtained analytically 

as 



fc- 1 fc 



and 



mf_ {k~l)K k fAKT^ 
2 2 \ TTa^ 



(17) 



ln(2). (18) 



We have checked, for the simplest case n — 2, that these 
expressions are very close to the results of the full (nu- 
merical) functional minimization of ^[p], for the relevant 
region of if > 20 and pi < lO^^cr^^. 

From Eqs. and we get f^™^ and w"-^ in Eq. 
(p^, for each of the local minima 



and 



(fc - 1) if 



fc In 



4^ 



1/2 



2nK^I'^ 



1/2 



ln(2) (19) 



(20) 



Both magnitudes become independent of the phase space 
unit volume P; they depend on the parameters of the 
interaction potential n, K and tr, the latter providing 
the only physical length scale in the model. 

We now use Eqs. (|l| and ^ in Eqs. (0 and |l|) 
to g et the MFA result for the total density p. In Fig. 
[I A the results of this MFA are compared with the exact 
results along the isotherm K — 2Q for n — 2 and 12. The 
qualitative trends are similar: the ideal gas is recovered 
at very low bulk density, but there is a smooth transition 
(around the CMC) to a regime in which /3fi has a slope of 
1/n, which corresponds to the ideal gas of full aggregates. 
In this regime the MFA goes nearly parallel to the exact 
curves, until the packing effects (neglected at this level of 
the MFA) become relevant. The CMC, still defined from 
nPn — pi, is shifted from the exact asymptotic result Eq. 

by 



l3{i^T^ - l^c) 



2n 



■ln(2) 



1 



2{n-l)' 



(21) 



In Fig. [[I A| we compare the the exact and the MFA re- 
sults for the CMC with n — 2 and 12. In the asymptotic 
regime the difference is independent of K. The MFA 
prediction for the CMC is shifted to larger values of the 
chemical potential: P{p^-^ — Pc) = 1/2, for n = 2, and 
In \/2 = 0.346 for large n. Within our approach this shift 
refiects the combination of the error made by the MFA 
in the estimation of (iE™^ (in the asymptotic regime Eq. 
( [TtI ) is 1/2 larger than the exact value), and the error 
associated with the MFA estimation of the cluster inter- 
nal entropy required to get O™^ and . In our simple 
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system we may compare the latter with the exact v„, as 
defined in Eq. (^) . In the Umit K ^ 1 we get 



en 
2K 



1/2 



(22) 



This exact result is independent of the optimal cluster 
size, n, and it is proportional to T^/"^ (through the pa- 
rameter K ^ T^^). The value given by the mean field 



approximation Eq. (12) shows the same dependence with 
the temperature and becomes identical to the exact value 
for n = 2. For larger n the result Eq. ( p2| ) goes like 2^/^ 
and saturates at large n with a discrepancy of with 
the exact result. As comparison, the identification of 
Vn with the thermal wavelength of the aggregate would 
lead to a qualitatively different dependence on n and T, 
Vn = K = h{2nl3/in m)y/^ ~ (n T)-^/"^, and with an 
unphysical dependence on the molecular mass and h. In 
our simple, one-dimensional, model the aggregates have 
a strong crystalline structure, described by the narrow 
Gaussian peaks in the MFA density distribution, which 
becomes very rigid a.t K 1. This character of the 
aggregates is responsible for the complete failure of t 
empirical recipe for u„ which was proposed by one of usi 
in terms of an aggregate compressibility. That empiri- 
cal recipe was used in a three-dimensional model, with 
much softer aggregates in which every molecule could 
move over the whole aggregate, like in a liquid drop. The 
application to our present model leads to a problematic 
formulation of u„ (as the square root of a negative num- 
ber) and should be fully discarded. 



tridimensional modelsH, and the consistent mean field re- 
sult for the density of clusters Eqs. (^ and ^5|) may be 
compared with the exact result. The comparison serves 
to validate our proposed approximation Eq. ( |l5|) for the 
prefactor Vn in the density of aggregates. As it should 
be expected, this prefactor is not related to the ther- 
mal wavelength of the aggregate, which would imply an 
unphysical dependence on the molecular mass and the 
wrong dependence on the temperature and the size of 
the aggregate. Our approximation reproduces the exact 
dependence on the temperature, and for large aggregates 
n 3> 1 it recovers (asymptotically) the exact indepen- 
dence on n. The numerical discrepancy between the ex- 
act and the MFA values of w„ should be regarded as in- 
herent to the mean field description of the aggregates. 
The thermodynamic consistency of the approximations 
for the exponent Qn and the prefactor w„ in the density 
of aggregates is the new feature of our proposal, and it 
could be used for any, more realistic, model of mesoscopic 
aggregates. 
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IV. CONCLUSIONS 



We have approached the study of ideal mixtures of 
molecular aggregates from a mean field approximation 
for the structure and internal free energy of the isolated 
aggregates. The general analysis in Section III A suggests 
a consistent approximation to get the missing factor in 
the relation Eq. (^ between the internal grand-potential 
energy of a cluster and its abundance in the equilibrium 
solution. The simple one-dimensional model of molecu- 
lar aggregation presented in Section II, provides a test 
for our approach. As in more realistic models, it includes 
hard-core repulsions and saturating attractive interac- 
tions, which set an optimal cluster size n. The exact so- 
lution of the model shows qualitative agreement with the 
behaviour of a system of amphiphilic molecules around 
the critical micellar concentration. Along an isotherm the 
system goes from an ideal gas of independent molecules 
(at very low chemical potential) to an ideal gas of op- 
timal clusters or micelles, at higher chemical potential. 
This transformation takes place over a range of fi which 
becomes narrower as n increases (although at a very 
slow rate because of the one-dimensional character of the 
model). The same model is then solved within a mean 
field approximation equivalent to what may be done in 
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